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^ Abstract. A procedure to obtain the eigenenergies and eigenfunctiones of a quantum 

spiked oscillator is presented. The originality of the method lies in an adequate 
' use of asymptotic expansions of Wronskians of algebraic solutions of the Schrodinger 

equation. The procedure is applied to three familiar examples of spiked oscillators 

Ph. 
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Spiked harmonic oscillators, i. e., harmonic oscillators to which a singular repulsion at 
the origin 



m ■ A 

, — , A > 0, a > 0, 

p . 

! has been added, have deserved a considerable interest since the publication of the 
^ ' pioneering paper of Klauder [1], three decades ago. A very extensive list of articles 
dealing with the topic can be found in more recent papers by Saad, Hall and Katatbeth 
[2] and by Liverts and Mandelzweig [3]. 
^ I Singular (at the origin) potentials were firstly considered in the context of collision 

I of particles the main interest being the adequate definition of the S matrix. Here we 
are concerned with the determination of the energy levels and the corresponding wave 
functions of a particle bound in a confining potential behaving at infinity like r", n > 2, 
and presenting a singularity of the type r~™, m > 2, at the origin. 

Besides of numerical integration of the Schrodinger equation by procedures adapted 
to the singular potential [5l[6l[7] and iterative methods of the Lanczos type |8], different 
approximate methods have been implemented to solve spiked oscillator problems. The 
usual WKB method was discarded from the very beginning in view of the results of 
Detwiler and Klauder ^ which also showed that conventional perturbative methods 
could not be used in the case of supersingular potentials. Harrell, in a very lucid 
paper [10], suggested a modified perturbation theory to a finite order applicable in 
that case. A perturbative study of nonsingular {a < 5/2) spiked harmonic oscillators 
in the two regimes of weak coupling (A ^ 1) and strong coupling (A ^ 1) has been 
carried out by Aguilera-Navarro et al [llj. They and other authors |12j have discussed 
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the connection between the expansions obtained in the two regimes. For the case of 
critical {a = 5/2) singularity, a large coupling perturbative expansion has been obtained 
|13] . One or more supersingular {a > 5/2) terms added to the harmonic oscilator 
potential have been treated also perturbatively by Guardiola and Ros [H] and by Hall 
and coworkers flbi [16] . Nonperturbative procedures have also been used: let us mention 
different implementations of variational methods [21 [131 [161 [T71 [HI [19] , matricial [20] and 
Hill-determinant [21] techniques, a conveniently modified WKB method [22], smooth 
transformations of solvable potentials [23], the pseudoperturbative shifted-/ expansion 
technique [23], and the more recently proposed quasilinearization method [3]. 

All these procedures provide with reasonably accurate values of the energies, at 
least for the ground state, but are not able to give the wave function in the vicinity of 
the origin. This fact, in addition to being unsatisfactory from an aesthetical point of 
view, may lead to incorrect expected values of certain operators. 

The first attempt to obtain the correct wave functions of spiked harmonic oscillators 
was done, to our knowledge, by Znojil [251 ES], who looked for solutions of the 
Schrodinger equation in the form of Laurent series multiplied by a non integer power of 
the variable. Substitution of such series in the differential equation leads to a recurrence 
relation satisfied by the coefficients; convergence of the series occurs only for a specific 
value of the non integer exponent. For any given energy, two independent solutions of 
that kind are obtained. They constitute a basis in the space of solutions of the differential 
equation for that value of the energy. The wave function, that can be written as a linear 
combination (with only one effective degree of freedom) of these two basic solutions, 
must be well behaved at the origin and at infinity. This double requirement can be 
satisfied, by adjusting the degree of freedom of the mentioned linear combination, only 
when the considered value of the energy is one of its eigenvalues. This is, schematically, 
the procedure followed by Znojil to find the energies and the wave functions of certain 
spiked harmonic oscillators. The procedure, rigorous in principle, presents the drawback 
that there is no possibility of evaluating the basic solutions neither at the origin nor 
at infinity. For this reason, Znojil obtained approximate values of the eigenenergies by 
requiring the vanishing of the wave function at a pair of points tq ^ 1 and r^o ^ ^■ 

The procedure that we present in this paper is much in the spirit of the Znojil's 
method, but differs from it in the manner in which the regularity of the wave function is 
imposed: our tq and are actually the origin and the infinity. Although the wave 
function cannot be calculated at these points, we are able to know its asymptotic 
behaviour in their vicinity and we may require that it be the correct one. The 
eigenenergies determined in this way are, therefore, exact save for errors inherent to 
the computation process. But these errors can be reduced by increasing the number of 
digits carried along the calculations. 

In the next Section we consider basic solutions of the Schrodinger equation for a 
very general three-dimensional spiked oscillator represented by the potential 




(1) 
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where the index q runs along a finite set of negative and positive integer and/or rational 
numbers, whose extremes are gj^jn < and gmax > 0. Moreover, in order to have a 
true spiked oscillator, we assume that both A{q-^[-^) and A(gniax) are positive. Section 
3 explains how to determine the behaviour at the origin and at infinity of a general 
solution. The requirement, in Section 4, of a regular behaviour at both singular points 
provides a quantization condition whose fulfilment determines the energy levels and the 
wavefunction. The procedure is applied to three popular cases of spiked oscillators in 
Section 5. Appendices A and B give details of the nontrivial steps of the method. 



2. Solutions of the Schrodinger equation 

The radial Schrodinger equation for the wave of angular momentum C 

£(£ + 1) 



with a potential given by (II]), can be easily written in the form 



R{r) = ER{r), (2) 



2 



d'^w 



^ ^^,+9iz)w = 0, (3) 

with a "potential" (including centrifugal and energy terms) that, multiplied by z"^, 
becomes a combination of positive and negative integer powers of z, 

2N 

9{z)= 9s z% M,N>0, (7_2M>0, g2N > 0. (4) 

s=~2M 

Obviously, z represents a power of the radial variable r, with exponent conveniently 
chosen, and w{z) is the reduced radial wave function R{r) multiplied by an adequate 
function of r. The origin and the infinity are the only singularities, of rank^ respectively 
M and A^, of the differential equation ([3]). We are interested in considering three pairs 
of independent solutions, namely: 

• Two Floquet or multiplicative solutions [281 EH], wi and W2, that, except for 
particular sets of values of the parameters in (jl]), have the form 



oo 



W,- = Z 



'3 



''3 

n=—oo 



E c„j2;", being ^ |c„jp < cx), j = l,2. (5) 



The indices Vj are not uniquely defined. They admit addition of any integer (with 
an adequate relabeling of the coefficients). In the general case, the indices Vj and 
the coefficients c„j may be complex. 

I We adopt the definition of ranks used in Ref. [37] , mainly because our work is based on the Naundorf 's 
treatment of the homogeneous linear second order differential equation. According to the definitions 
used in Ref. [3S], the singularities should be of ranks M + 1 and + 1. 
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• Two Thom^ formal solutions, W3 and w^, that have the nature of asymptotic 
expansions for 2; — )■ 00, 

^ ^ Um. ^ , ao,k ^0, = 3, 4. (6) 

p=l P J m=0 

It is usual to say that these two expansions are associated to each other. 

• Two r/iome formal solutions, and we, asymptotic expansions for z — )■ 0, of the 
form 

(Mo \ 00 

E — ^"' P'' E^™,^-^"' V^^O, / = 5,6. (7) 
q=l 1 J m=0 

Also these expansions are associated. 

The determination of the indices Vj and the coefficients c„j of the multiplicative 
solutions is rather laborious. By substitution of ([5]) in (|3]) one obtains the infinite set of 
homogeneous equations for the coefficients 

2N 

(n + z/j)(n-l + z/j) c„j - E fl'sC„-sj = 0, n = . . . , -1, 0, 1, . . . , (8) 

s=-2M 

that can be interpreted as a nonlinear eigenvalue problem, where the eigenvalue v must 
be such that 

00 

E |Cn,jf <00. (9) 
n=— 00 

In Appendix A we recall the Newton iterative method to solve that problem. In general, 
two indices, z/i and z/2, and two corresponding sets of coefficients, {c„,i} and {c„,2} are 
obtained, but for certain sets of values of the parameters Qs only one multiplicative 
solution appears. Any other independent solution must include powers of the variable 
multiplied by its logarithm. Such logarithmic solutions cannot correspond, usually, 
to the practical system that one tries to describe and are, therefore, to be discarded. 
We will assume, from now on, that the parameters Qs are such that ([3]) admits two 
independent multiplicative solutions. 

In what concerns the formal solutions W3 and W4, the exponents a and /i and the 
coefficients am must be such that the expansions in the right hand side of (E]) satisfy the 
differential equation. One obtains in this way 



N \ ^ iV 



E "p^^ + E(^'^i)"p^^ - 9{z) 



00 I N \ 00 

am Z'"^ + 2 E "P^^ E (-"^ + /^) Z~ 
m=0 \,P=1 / rn=0 



+ Y,{-m+n){-m-l+n)amZ~'"'' = Q. (10) 



m=0 



§ Although the credit of these solutions is attributed to Thome, they were used before by Fabry and 
by Poincare. We are grateful to Prof. Alberto Griinbaum for illustrating us about this fact. 
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Cancellation of the powers z"^^ , z"^^'^, . . . , z^~^^ inside the bracket of the first term in 
f irUj) produces a system of equations 

(ctAr)^ - g2N = 0, 

2 aN 0!N~1 - g2N-l = 0, 

2aAraAr_2 + (aAT-l)^ - 5'2Ar-2 = 0, (11) 



2 Oat «! + ... - g^+i = 0, 

which can be solved successively. There are two sets of solutions, {op.a} and {ap,4}, that 
obviously verify 

ap,3 = -o;p,4, p=l,...,N. (12) 
Let us assign the labels 3 and 4 in such a way that 

«iV,3 = -^/g2N, OiN,i = +^/g2N■ (13) 

Accordingly, the formal solution w^, presents, from the physical point of view, the 
adequate behavior at infinity, on the positive real semiaxis, whereas W4 should be 
rejected. Cancellation of the coefficient of z^ in the left hand side of (fTOl) gives the 
two exponents fi^ in terms of the previously obtained ap^k- Notice that 

^l3 + ^^4 = -N + 1. 

Finally, cancellation of the coefficient of z^, n = N — 1, N — 2, . . . , — 00, implies a 
recurrence relation for the coefficients am,k that allows one to obtain all of them starting 
with an arbitrarily chosen ao.fc 7^ 0. 

Analogously, by requiring the expansions in the right hand side of (I7j) to satisfy the 
differential equation, one obtains for the formal solutions 1^5 and Wq 

^ A/ ] 00 / M \ 





m=0 \<7=1 / m=0 



+ £ (m+p)(m-l+p) h^z^ = 0. (14) 



m=0 

'2M ^-2M+l ■ 



Now, cancellation of the powers z , z ^^^+1^ . . . ,z inside the bracket of the first 
term in f|T^ gives the system of equations 

ihlY - 9-2M = 0, 

2 Pm Pm-I — g-2M+l = 0, 

2 Pm (^M-2 + {(^M^lf - 9-2M+2 = 0, (15) 



2(3MPi + ...-g-M-i = 0. 
The two sets of solutions, {/3g,5} and {f3gfi}, verify 
/9,,5 = -/3,,6, g=l,...,M. 



(16) 
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If the labels 5 and 6 are assigned in such a way that 

/3m,5 = —■\/g-2M, /3m,6 = +y/g-2M, (17) 

the formal solution vanishes at the origin, whereas wq diverges. Cancellation of the 
coefficient of in the left hand side of (fl^ allows one to obtain the two exponents 
pi in terms of the previously calculated (3g^i. They obey the relation 

P5 + P6 = M + 

The coefficients bm,i can be obtained, starting with an arbitrarily chosen Bq i, by making 
use of the recurrence relation steming from the cancellation of the coefficient of z", 
n = -M + 1, -M + 2, . . . , +00, in the left hand side of (lT4l). 

3. The connection factors 

Any solution w of the differential equation ([3]) can be written as a linear combination 
of the two multiplicative solutions, 

W = C1W1 + (2 W2. (18) 

Its behavior in the neighbourhood of the singular points can be immediately written 
if, besides the coefficients Ci and C2, one knows the behaviour of the multiplicative 
solutions, that is, if one knows the connection factors T of their asymptotic expansions, 

Wj ~ Tj_3 W3 + Tj 4 W74, for 2; — )• 00, j = 1, 2, (19) 
Wj ~ T,- 5 W5 + Tjfi We, for z 0, j = 1, 2. (20) 

These connection factors are, obviously, numerical constants, but their values depend on 
the sector of the complex plane where z lies. This fact, known as "Stokes phenomenon" 
|30] . introduces a slight complication. As it is well known, the connection factor 
multiplying any one of the asymptotic expansions in the right hand sides of ( |T9l) and 
( 120|) takes different values in the sectors of the complex 2;-plane separated by a Stokes 
ray of the associated expansion. On the ray, the value of the connection factor is the 
average of those two different ones. We are interested in the behaviour of w{z) in the 
vicinity of the origin and at infinity on the ray arg^; = 0. Since this is a Stokes ray for 
the expansions and wq, the values of Tj^a and T^-^s for arg2; = are the respective 
averages of their values for z just above and below the positive real semiaxis. 

The problem of finding the connection factors T has been considered by several 
authors. Most of them |3T] refer to a differential equation for which the origin is an 
ordinary or a regular singular point and the infinity is an irregular singular one. As 
far as we know, only the procedure developed by Naundorf [27] becomes applicable also 
to the present case of both the origin and the infinity being irregular singular points. 
In a former paper [32] we have suggested a modification of the Naundorf 's procedure, 
improving notably its performance, and applied it to find the bound states of anharmonic 
oscillators, whose Schrodinger equations present an irregular singular point at infinity, 
the origin being an ordinary or a regular singular one. Some mathematical questions 
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left aside in that paper have been tackled in a posterior one [53]. As we are going to 
show, our improvement of the Naundorf 's method is easily applicable in the case of the 
origin being also an irregular singular point, as it happens in the Schrodinger equations 
of spiked oscillators. 

Our procedure to calculate the connection factors makes use of the fact that they 
can be written as quotients of Wronskians of the solutions presented in section 2. We 
adopt the usual definition of Wronskian of two functions f{z) and g{z), namely 

We benefit from the fact that equation ([3]) is a second order linear differential one where 
the first derivative term is absent and, therefore, the Wronskian of any two of their 
solutions is a constant. Then, it is immediate to obtain 

The denominators of the quotients in the right hand sides of these equations can be 
computed trivially. One obtains 

W[w3, W4] = - W[w4, W3] = - 2aN,3 ao,3 ao,4, (23) 
WK, We] = - W[we, W5] = 2/3m,5 &o,5 &o,6- (24) 

The calculation of the numerators is not so easy. Direct computation of any one of them 
gives a certain power of z times an infinite series of negative and positive powers of z. 
This does not seem to be an adequate expression to determine the (independent of z) 
value of the Wronskian. Instead, we have devised a trick that allows one to obtain each 
one of those numerators. Appendix B contains a detailed description of the procedure. 

4. The quantization condition 

We have already mentioned that any solution of (E]) can be written as a linear 
combination of the two multiplicative solutions, like in f|T8|) . According to f l20|) . its 
behaviour in the neighbourhood of the origin is given by 

w{z) ~ (Ci Ti,5 + C2 T2,5)w5{z) + (Ci Ti,6 + C2 T2,e)we{z) for z 0. (25) 
Since Wq diverges as z goes to zero, the coefficients Ci and (2 must be such that 

Tl,6Cl+T2,6C2 = 0. (26) 

On the other hand, at infinity, 

w{z) ~ (Ci Ti,3 + C2 T2,3)w3{z) + (Ci Ti,4 + (2 ^2,4)^4(2) for z-)- 00. (27) 
The correct asjTiiptotic behaviour occurs when 

Ti,4 + ^2,4(2 = 0. (28) 
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The fulfilment of both requirements fl26p and fl28l) implies the quantization condition 

Tl,6T2,4 - Ti,4T2,6 = 0. (29) 

For given potential ([T]) and angular momentum, the connection factors are functions 
only of the energy, through one of the parameters Qg in (jlj). Therefore, the left hand 
side of (!29|) is but a function of the energy whose zeros are the eigenenergies of the 
spiked oscillator. 

The problem of finding the wave functions is now immediately solved. Any 
nontrivial solution {Ci,C2} of the system of equations ( l26i) and ( l28l) gives the 
(unnormalized) physical solution 

^i^phys(2;) = Cl + C2W2{z). (30) 

For the computation of Wphys in the neighbourhood of the origin or for large values of 
z, the asymptotic expansions 

Wphysiz) ~ (Cl Ti,5 + (2 T2,5)t(75(z), Z ^ 0, (31) 

^i^phys(2;) ~ (Cl Ti-i + 4 T2,3)iy3(z), Z OO, (32) 

should be used instead of ( |30l) . that suffers from strong cancellations in the right hand 
side for z in those regions. So, there is no difficulty to normalize Wphys properly. 



5. Some examples 

In order to illustrate how our procedure may serve to treat spiked oscillators, we are 
going to apply it to a few cases that have already been considered, by having recourse 
to different approximations, by other authors. The numerical results given below have 
been obtained by using double precision FORTRAN codes. In our opinion, the digits 
shown in the quoted values of the energies are correct. Our procedure could provide 
additional correct digits but it would require a higher precision arithmetic. This need is 
more evident in the case of high values of the potential parameters, angular momentum 
different from zero, or excited states. 

Following a common practice and to facilitate comparison of our results with those 
of other authors, we assume that r is a dimensionless variable that represents a distance 
in a given scale. Analogously, we use dimensionless symbols, E and Ag to represent 
the energy and the intensities of the different terms of the potential, that are assumed 
expressed in adequate units. 

5.1. Potential V{r) = A2 + A_4 

The first example to be considered is a three dimensional spiked harmonic oscillator, of 
potential 

y(r) = A2r2 + A_4r-^ (33) 
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which has been most discussed by other authors. We will assume, without loss of 
generality, that 

^2 = 1. 

The radial Schrodinger equation ([2]), written in terms of the variable 

z = r , 

turns into Eq. ([3]) for the function 

w{z) = R{r), 

with 

g{z) = A_4 + C{C + l)-Ez'^ + z\ 

where, obviously, C represents the angular momentum quantum number. So, the ranks 
of the singularities at the origin and at infinity are, respectively, 

M = l, N = 2. 

The coefficients c„ j- of the Floquet solutions obey the recurrence relation (omitting the 
second subindex j) 

The Thome solutions at infinity have exponents 

-1 + E -1-E 

a2,3 = -a2,4 = -1, ai,3 = -ai,4 = 0, /is = fi4, = , 

and coefficients a^^k {k = 3,4) obeying the recurrence relation (omitting the second 
subindex k) 

2a2mam = [{m — 2 — /i)(m — 1 — yu) — £(£ + 1)] 0^-2 — ^-4 dm-i- 
For the Thome solutions at the origin the exponents are 

/3l,5 = -/3l,6 = -^-^ , P5= P6 = 1, 

and the coefficients bm,k, ik = 5, 6) satisfy the recurrence relation (omitting the second 
subindex k) 

2f3imbm = [m{m - 1) - £(£ + 1)] bm-i + E h^^^ - h^-b- 

The definitions of Floquet and Thome solutions given above fix them up to arbitrary 
multiplicative constants. To avoid ambiguities in the definition of the connection factors, 
we assume from now on that those arbitrary constants have been chosen in such a way 
that 

Co,l = Co,2 = 1, ^0,3 = 00,4 = 1, &0,5 = ^0,6 = 1- 

Then, according to Eqs. ([23]) and (^^, we have for the denominators in Eqs. (1^ and 

WK, H = - W[^i;4, ws] = 2, (34) 
W[w,, We] = - W[we, w,] = -A^l- (35) 
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Our procedure gives for the numerators in the same equations (for j = 1,2) 




(36) 
(37) 
(38) 
(39) 



where we have abbreviated 



5r^ = (i/,+/ifc + l)/2 



= 3,4, 



(40) 



oo 



-m+s— l,j 



{m + 2s + l + Uj-fik)Cm+s+l,j) , 



A; = 3,4, 



(41) 



s=0 



oo 



■m—s- 



:+2j + (2s - m - Uj) c. 



■m—s- 




1 = 5,6. 



(42) 



s=0 



We report in Table 1 our ground-state energies of the spiked oscillator of potential 
( 133|) . for several values of the parameter A_4. The quoted energies should be compared 
with those in Table 3 of Ref . [H] , that collects results obtained by several authors with 
different methods. Obviously, our procedure is superior to variational ones, as it provides 
with four more correct digits in the energies for the values of y4_4 considered. In fact, 
our double precision FORTRAN calculations give results with an accuracy comparable 
to the excellent one reached by Buendia et al. [7] thanks to an analytic continuation 
method or by Roy [6] by using a generalized pseudospectral method. (Notice that, due 
to a different definition of the Schrodinger operator, our energies should be divided by 2 
if they are to be compared with those of Table 4 of Ref. [6].) As mentioned above, larger 
values of A^^ require a more precise arithmetic to the same accuracy in the results. For 
completeness, we give also in the same table the indices Uj of the Floquet solutions. We 
have already mentioned the existing ambiguity in the values of the indices. To avoid 
it, we assume that the indices vary continuously with the parameter y4_4 and fix their 
integer part in such a way that i^i = (and, consequently, 1^2 = 1) for A_4 = 0, as 
it should be for a particle of zero angular momentum. As increases, the index ui 
increases also until it reaches the value 0.5 for A_4 = 0.1305... and an eigenenergy 
E = 3.6454.... Notice that, for those critical values of the parameters, ui = 1/2 and 
only one Floquet solution of the form (|5]) is obtained. Any other independent solution of 
the differential equation must contain logarithmic terms. Our procedure, in its present 
form, is not applicable in this case. A different implementation of the basic idea would 
be necessary. If A_4^ continues increasing above the critical value, the real part of ui 
remains equal to 0.5, whereas its imaginary part increases. Obviously, the two Floquet 
solutions are then complex conjugate to each other. This fact allows one to simplify the 
calculations in our procedure. 

In order to illustrate the trend of the coefficients of the Floquet and Thome 
solutions, we show, in Tables 2 to 7, the most relevant of them for two particular 
cases. Tables 2, 3 and 4 correspond to the ground state of the example discussed by 
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Table 1. Ground-state energy of a spiked oscillator of potential ([33l) for several values 
of the intensity A^^ of the singular repulsive potential. The index vi of one of the 
Floquet solutions, wi, is also given. The index 1/2 of the other Floquet solution, W2, is 
immediately obtained from the relation 1^2 — ^ — i^i- 



A_4 


E 






0.0001 


3.0222745087 


0.201485000573E-03 


0.001 


3.0687631709 


0.204586237797E-02 


0.005 


3.1483523083 


0.104967473634E-01 


0.01 


3.2050674951 


0.213850813448E-01 


0.1 


3.5755519912 


0.270240464647E+00 


0.4 


4.0319714400 


0.5 - 


i0.606083134346E+00 


1 


4.4941779834 


0.5 - 


f i0.950268234562E+00 


10 


6.6066225120 


0.5 - 


t- i0.203793867918E+01 


100 


11.265080432 


0.5 - 


f- i0.412681646514E+01 



Znojil [25], namely that of A_4 = 0.4. In this case, the energy and the indices of the 
Floquet solutions are 

E = 4.0319714400 , z/i = = 0.5 + z 0.606083134346, 

and the connection factors 



Ti 3 = T2 3 = 0.363722440835 + i 1.049440656062 , 



Ti 4 = T2 4 = - 0.378572197756 + i 0.137728550255 



Ti 5 = T2 5 = 0.520935174155 + i 1.431885933657 , 



Ti,6 = T2,6 = 0.436272922113 - i 0.158720681109 . 

The coefficients of the Floquet solutions can be seen in Table 2, and those of the 
Thome solutions at infinity and at the origin in Tables 3 and 4, respectively. The 
second particular case, illustrated in Tables 5, 6 and 7, is that of A_4 = 1 and angular 
momentum C = 2, considered by Aguiler a- Navarro and Ley Koo Table II]. For the 
energy and the indices of the Floquet solutions we have obtained 

E = 7.2242871639 , i^i = -2.083592228877 , U2 = 3.083592228877 , 

and for the connection factors 

Ti,3 = - 0.7041938314 , T2,3 = 0.9280831701 , 

Ti,4 = - 2.4617003408 , T2,4 = - 0.0462330978 , 

Ti^5 = - 7.7185228435 , T2,5 = 0.1449612757 , 

Ti,6 = 7.9925473468 , T2,6 = 0.1501077191 . 

The ambiguity in the integer part of the indices of the Floquet solutions has been 
resolved in this second particular case bearing in mind that, since £ = 2, it is desirable 
that ui — )■ —2 and 1^2 ^ 3 for — )■ 0. 

Once the connection factors are known, the wave function is easily obtained 
following the procedure indicated in Section 4. 
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Table 2. Some coefficients of the Floquet solution wi, normalized in such a way that 
co^i = 1, for the particular case of A_4 = 0.4, £ = and energy E = 4.0319714400. 
The coefficients with odd first label are equal to zero. The other Floquet solution, W2, 
has coefficients Cm,2 complex conjugate of the Cm,i- 



-10 -0.142317651396E-23 

-9 -0.100263927421E-20 

-8 -0.429957466419E-18 

-7 0.167970940662E-17 

-6 0.156002041860E-12 

-5 0.120892275974E-09 

-4 0.492753885758E-07 

-3 0.115445648553E-04 

-2 0.144400180631E-02 

-1 0.780262828557E-01 

1 

1 -0.760452558754E+00 

2 0.199603204168E+00 

3 -0.348689735973E-01 

4 0.433386333626E-02 

5 -0.431181023059E-03 

6 0.349885697033E-04 

7 -0.243278698241E-05 

8 0.146155963833E-06 

9 -0.779445127491E-08 
10 ().37()(ir)48(i99.^)7E-()9 



0.690261395778E-23 
0.697827393080E-20 
0.569639303294E-17 
0.365781757971E-14 
0.178667148961E-11 
0.634808135173E-09 
0.154076769690E-06 
0.232262850896E-04 
0.184616207302E-02 
0.537272549009E-01 

0.536951488093E+00 
-0.202477697179E+00 

0.453398033971E-01 
-0.673686840525E-02 

0.781829457162E-03 
-0.725041819983E-04 

0.570826186141E-05 
-0.385089511701E-06 

0.229774575047E-07 
-().;i2LS(i21,S4787E-08 



Table 3. The first coefficients of the two Thome solutions at infinity for the particular 

case of = 0.4, C = and energy E = 4.0319714400. We have skipped over the 
coefficients with odd first label, since they are equal to zero. 



m 


a2m,3 


a2m,4 





1 


1 


1 


-0.195556745811E+00 


0.221154246581E+01 


2 


0.675581627248E-01 


0.683622152000E+01 


3 


-0.552411742467E-01 


0.278260729988E+02 


4 


0.865891090477E-01 


0.140764673339E+03 


5 


-0.212297950675E+00 


0.851787239516E+03 


6 


0.709506391974E+00 


0.6001543154448E+04 


7 


-0.305385835605E+01 


0.482637790136E+05 


8 


0.160735596654E+02 


0.436250311536E+06 


9 


-0.100168048970E+03 


0.437841970329E+07 


10 


0.721888409401E+03 


0.483138693754E+08 
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Table 4. The first coefficients of tlie two Thome solutions at the origin for the 
particular case of A_4 = 0.4, £ = and energy E = 4.0319714400. 



71 i 


''IN ■ ) 







1 


1 


-1 
1 


U 


U 


2 








3 


-0.106251776760E+01 


0.106251776760E+01 


4 


0.251998215001E+01 


0.251998215001E+01 


5 


-0.781076937396E+01 


0.781076937396E+01 


6 


0.314392488782E+02 


0.314392488782E+02 


7 


-0.150276962631E+03 


0.150276962631E+03 


8 


0.834637749300E+03 


0.834637749300E+03 


9 


-0.528962618019E+04 


0.528962618019E+04 


10 


0.376836341629E+05 


0.376836341629E+05 



Table 5. Coefficients of the Floquct solutions wi and normalized in such a way that 
co,i = Co, 2 = 1, for the particular case of A_4 = 1, £ = 2 and energy E = 7.2242871639. 
The coefficients with odd first label are equal to zero. 



C2n,l 



C2n,2 



-10 0.418079159329E-21 

-9 0.210620217409E-18 

-8 0.879231161250E-16 

-7 0.298162410599E-13 

-6 0.801421283847E-11 

-5 0.165459858432E-08 

-4 0.251717006169E-06 

-3 0.266341461331E-04 

-2 0.179770550637E-02 

-1 0.668756857204E-01 

1 

1 0.906280257537E+00 

2 0.180182069867E+01 

3 -0.209931807620E+01 

4 0.728418847644E+00 

5 -0.150546747680E+00 

6 0.220037264283E-01 

7 -0.249248011094E-02 

8 0.230170752863E-03 

9 -0.179507697993E-04 
10 0.121105001156E-05 



0.508287996884E-16 
0.151012739304E-13 
0.349503786349E-11 
0.607371865088E-09 
0.753905871024E-07 
0.621798117715E-05 
0.303689920293E-03 
0.705627796961E-02 
0.404457498845E-01 
-0.120970341064E+00 
1 

-0.489421426933E+00 

0.121731165609E+00 
-0.202632214005E-01 

0.253270584477E-02 
-0.253385266560E-03 

0.211291907657E-04 
-0.151036531706E-05 

0.944714256302E-07 
-0.525254239416E-08 

0.262829669113E-09 
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Table 6. The first coefRcients of the two Thome solutions at infinity for the particular 
case of A_4 = 1, £ = 2 and energy E = 7.2242871639. We have skipped over the 
coefficients with odd first label, since they are equal to zero. 



m 


0'2m,3 









1 




1 


1 


-0.143323523200E+00 


0.375546710516E+01 


2 


0.197417671158E- 


-01 


0.174648768948E+02 


3 


-0.483024944799E- 


-02 


0.985368170788E+02 


4 


0.281201497378E- 


-02 


0.653978688776E+03 


5 


-0.344425647424E- 


-02 


0.499199894361E+04 


6 


0.705309939217E- 


-02 


0.430838076974E+05 


7 


-0.207487524900E- 


-01 


0.414831591842E+06 


8 


0.802543608266E- 


-01 


0.440833793663E+07 


9 


-0.386209034478E+00 


0.512489320365E+08 


10 


0.22278847;Ui()lE4 


-01 


O.G4(i9,S319,sr,G,SE+09 



Table 7. The first coefRcients of the two Thome solutions at the origin for the 

particular cast- of A_i = 1. C = 2 and cuvrgv E = 7.2242871639. 



m 


&m,5 







1 


1 


1 


3 


-3 


2 


3 


3 


3 


-0.120404786066E+01 


0.120404786066E+01 


4 


-0.180607179098E+01 


-0.180607179098E+01 


5 


0.461214358197E+00 


-0.461214358197E+00 


6 


0.524369089808E-01 


0.524369089808E-01 


7 


0.101141803810E+01 


-0.101141803810E+01 


8 


-0.344418092084E+01 


-0.344418092084E+01 


9 


0.125072805387E+02 


-0.125072805387E+02 


10 


-0.528728562620E+02 


-0.528728562620E+02 



5.2. Potential V{r) ^ A2r^ + r'"^ + ^-6 r'^ 

Our second example is a three dimensional spiked harmonic oscillator which has received 
also considerable attention and which becomes quasi-exactly solvable for certain sets of 
parameters. The potential is 

y(r) = A2r' + A_4r-^ + A-6r-^ (43) 

Following a common practice, we will assume 

^2 = 1. 

We can reduce the rank of the singularities of the radial Schrodinger equation by using 
a variable 
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Equation ([2]) turns in this way into Eq. ^ for the function 

w{z) = r^/^R{r), 

with 

g{z) = i [A^6 z-^ + A_4 z-^ + CiC + l)-3/A-Ez + z^) . 
Now the ranks of the singularities at the origin and at infinity are, respectively, 
M = l, N = l. 

This fact reveals that we are dealing with a double confluent Heun Equation [31]. In a 
recent paper [55], we have detailed the algorithm resulting from the application of our 
method to an equation of this type. For the coefficients c„j- of the Floquet solutions we 
have the recurrence relation (omitting the second subindex j) 

A_G c„+4 + A_4 c„+3 + + 1) - 3/4 - 4(n + 2 + i^)(n + 1 + z/)] c„+2 - E c„+i + c„ = 0. 
The Thome solutions at infinity have exponents 

"1,3 = -"1,4 = -1/2, fl3=E/4:, iJ,^=-E/4:, 

and the recurrence relation for the coefficients am,k {k = 3,4) is (omitting the second 
subindex k) 

8 aim am = [4(m - fi){m - 1 - fi) - C{C + 1) + 3/4] a„_i - A_4a„_2 - ^-e^m-s- 

For the Thome solutions at the origin, the exponents are in this second example 

/3i,5 = -/3i,6 = -A'-_I^J2 , P5 = 1 + A.,/AA'l^, , p^ = l- A.JAA^I^ , 

and the coefficients hm,i {I = 5,6) obey the recurrence relation (omitting the second 
subindex /) 

8 /3i m 6^ = [4(m - 1 + p)(m - 2 + p) - £(£ + 1) + 3/4] hm-i + E hm-2 - b^-s- 
By choosing, as in the first example, 

Co,l = Co,2 = 1; O'0,3 = ^0,4 = 1; ^0,5 = ^0,6 = 1- 

we have, from Eqs. fl23l) and flMl) . 

W[w3,w4 = -W{w,,W3] = l, (44) 
W[w,, We] = - W[we, w,] = -A^l- (45) 

For the numerators in Eqs. (l2Til and ( |22l) our procedure gives in this second example 
(for J = 1,2) 

WK, ti;3] = 2"+'^'"'' r(n + 1 + ^(^'3)) 7(^'3), (46) 
«;4] = (-1)" cos(7r5(^'^)) 2"+^*''*' r(n + 1 + (J^^-^)) 7(^'^\ (47) 
WK, ^5] = (Ai_4V2)-"-^'^-^'' r(n + 1 + 5(^-'5)) 7^;f \ (48) 
WK-,W6] = (-1)" cos(7r5(^'6)) (^^5/2)"""'^'''" r(n + 1 + 5(^'^))7i'n\ (49) 
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Table 8. Ground-state energy of a spiked oscillator of potential (|43l) for several values 
of the intensities A^q and A_4. The column headed by E contains the results obtained 
with our method. For comparison, we show, in the column headed by Em, three set 
of values taken from the literature, namely from Refs. [7], [6] and [2], respectively. 
(Notice that, due to a different definition of the Hamiltonian, a factor 2 has been 
applied to the results of Ref. [6] and, consequently, the last digit may oscillate by one 
unit.) The last column shows the indices vi and 1^2 of the Floquet solutions. 



A_6 




E 


Eiit 




0.001 





3.27985582592 


3.279855825921856 


0.249216175554 


0.0025 





3.35391931711 


3.353919317108725 


0.247958538878 


0.01 





3.50545227600 


3.505452275995097 


0.241137578178 


1 





4.65993996957 


4.659939969573538 


i 0.466911061788 


10 





6.00320902890 


6.003209028895745 


i 0.895534935089 


0.005 





3.42288418426 


3.42288418426 


0.245761020193 


0.05 





3.76554020606 


3.76554020604 


0.198535942381 


0.5 





4.38790906027 


4.38790906026 


i 0.337261268644 


5 





5.51315901419 


5.51315901418 


i 0.768433078693 


1 


10 


6.67905366445 


6.679054 


0.5-^ 1.00539309301 


10 


1 


6.14012287178 


6.140123 


i 0.896525791611 


10 


10 


7.13826093998 


7.138261 


0.5 -i 0.320864634688 



with the abbreviations 

50'^) = uj+iJk, k = 3, 4, = -Uj -pi, 1 = 5,6, (50) 

00 

lm''^ = J2^^,k{ai,kCm+s,j-{m + 2s + l + Uj-fXk)Crn+s+l,j), = 3, 4, (51) 

s=0 

00 

7m'^ = ^bs,l (- (3i^i Cm-s+2,j + {2s -m-1 - Uj + pi) Cm-s+l,j) , / = 5, 6. (52) 

The connection factors are then obtained immediately from Eqs. f l2T|) and fl22|) . 

The ground state energies of the spiked oscillator of potential f H3|l obtained with 
our method, for several values of A^q and A_4, are shown in Table 8. The chosen 
values of the parameters allow to compare our results with the extremely precise ones of 
Buendia et al. [7] and of Roy [6]. The concordance is remarkable. We have considered 
also values of A_q and v4_4, both different from zero, allowing comparison with the 
variational results of Saad, Hall and Katatbeth |2]. In the same table we give also the 
indices vi and z/2 of the Floquet solutions. As in the preceding example, the physical 
solution Wphys can be obtained easily following the steps indicated in Section 4. Going 
back to the original variable r and reduced radial wave function R{r) is trivial. 

5.3. Potential V{r) = + Ar'^/^ 



As a third example, we have chosen the spiked oscillator whose potential 

V{r) =r^ + Xr-^/^ (53) 
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presents a critical singularity. The Schrodinger equation ([2]) adopts the form for the 
variables 

z = r^/^ and w{z) = r~^^^R{r), 
g{z) being now 

g{z) = 16 (a z-^ + C{C + 1) + 15/64 - E z^ + z^^ 
The singularities at the origin and at infinity have ranks 

M = l, A^ = 8. 

The recurrence relation for the coefficients c„j of the Floquet solutions ([5]) is (omitting 
the subindex j = 1, 2) 

Ac„+i8 + [C{C + 1) + 15/64 - (n + 16 + z/)(n + 15 + z/)/16] Cn+w - Ec^+s + = 0. 
For the exponents of the Thome solutions at infinity we have 

«8,3 = — ^8,4 = —4 , 0:7 J = a6,i = • • • = aij = , H3 = —7/2 + 2E, = —7/2 — 2E, 

and for their coefficients am,k (omitting the subindex A; = 3, 4) 

2 ag m Om, = [(m — 8 — /i)(m — 7 — /i) — 16 £(£ + 1) — 15/4] ~ ISA cim-io- 

The Thome solutions at the origin have exponents 

/3l,5 = -A,6 = -4Al/^ P5 = pe = l, 

and coefficients bm.i that satisfy the recurrence relation (omitting the second subindex 
/ = 5,6) 

2(3imbm = [m{m - 1) - 16C{C+1) - 15/4] + lQEbm-9 - lQb^~i7- 
If we choose, as in the preceding examples, 

Co,l = Co,2 = 1, ^0,3 = CtO,4 = 1, ^0,5 = ^0,6 = 1; 

we have for the denominators in Eqs. (!2T|) and (!22|) 

>V[W3,H = ->V[W4,W^3]=8, (54) 

>V[«;5, We] = - W[we, w,] = -8 \'/\ (55) 
The numerators in the same equations (for j = 1,2) are given by 

WK-, ws] = E 2"+^"''' r(n + 1 + 5i^'^)) 7^1, (56) 

L=0 

WK-, u;4] = (-1)" E cos(vr5i^'^)) 2"+^i^''' r(n + 1 + 6^^) ^t^, (57) 

L=0 

WK-, H = (4Ai/2)-"+'^^+i r(n - z.,) 7^;f\ (58) 
WK-,w;6] = (-1)"-' cos(7rz/,)(4Ai/2)-"+^^+ir(n-i/,)7^;f\ (59) 
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Table 9. Ground-state energy of the critical spiked oscillator of potential (|53|) for 
several values of the intensity A. The column headed by E shows our results. For 
comparison, we show, in the column headed by Em, values taken from the literature, 
namely from Refs. [7] (superscript a) and [13 (superscript b), obtained by the analytic 
continuation method and by numerical integration, respectively. The last column shows 
the index i^i of one of the Floquet solutions. For the index of the other one we have 



taken i>2 












A 


E 


Em 








0.001 


3.00401125101 


3.004011251013044^ 


0.5 H 


hi0.244567376746E- 


-04 


0.005 


3.01914010728 


3.019140107276879^ 


0.5 H 


hi0.612956032070E- 


-03 


0.01 


3.03672947263 


3.036729 b 


0.5 H 


hi0.245895012676E- 


-02 


0.05 


3.15242944140 


3. 152429 b 


0.5 H 


hi0.625334357268E- 


-01 


0.1 


3.26687302611 


3.266873026113018^ 


0.5 H 


hi0.248671350579E4 


-00 


0.5 


3.84855317229 


3.848553 b 


0.5 H 


hi0.196189243685E4 


-01 


1 


4.31731168925 


4.317311689247366^ 


0.5 H 


hi0.288463702918E4 


-01 


2 


4.98613573609 




0.5 H 


hi0.389836484474E4 


-01 


5 


6.29647263890 


6.296472 b 


0.5 H 


hi0.599616514206E4 


-01 


10 


7.73511110349 


7.735111103489141^ 


0.5 H 


hi0.813625698416E4 


-01 


20 


9.70940409621 




0.5 H 


hi0.110399234356E4 


-02 



where we have abbreviated 

Sl'^ = i^j +f^k + L)/8 , A; = 3, 4, (60) 

oo 

Im''^ = XI "^.fc {0is,kCm+s-7,j - {m + 2s + 1 + Uj - fik) Cm+s+l,j) , k = 3,4, (61) 

s=0 

Im^^ =J2^s>^('~ (^^>l'^rn-s+2,j + i'2s-m-Uj)Cm-s+l,j), 1 = 5,6. (62) 

s=0 

The ground state energies obtained by using our procedure are reported in Table 
9. We have taken for A several values already considered by other authors, whose 
results are shown also in Table 9 in order to facilitate comparison. As it can be 
seen, our result are considerably more accurate than those obtained by numerical 
integration of the Schrodinger equation [13], although we do not reach, with our double 
precision FORTRAN codes, the impressive accuracy of the results of Buendia et al. 
[7] obtained with the analytic continuation method. Our procedure, however, can 
provide with results with as many correct digits as desired if a sufficiently precise 
arithmetic is used. And, very important, the wave function is obtained in a very 
convenient form (asymptotic expansions and Laurent series) for algebraic manipulations 
like normalization or computation of expected values. 
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Appendix A. Indices and coefficients of the Floquet solutions 

In Sec. 2 we have referred to the Floquet or multiplicative solutions ([5]) of Eq. (3) 
and mentioned that their indices Uj and coefficients c„j must obey the infinite set of 
homogeneous equations ([8]) and the condition (Q. (Along this Appendix we will omit, 
for brevity, the subindex j in Uj and c„j.) As it has been already said, we face a nonlinear 
eigenvalue problem. Algorithms to solve finite order problems of this kind have been 
discussed by Ruhe [36]. Obviously, the condition requires 

lim |c„| = 0. (A.l) 

n— >-±oo 

This allows us to truncate our infinite problem (jH]) by restricting the label n to the 
interval —Ai < n < A/", both Ai and JV being positive integers large enough to guarantee 
that the solution of the truncated problem does not deviate significantly from that of 
the original infinite one. Then, the Newton iteration method can be applied. It consists 
in moving from an approximate solution, {u^'^\c^^}, to another one, c^^'^^'^}, by 

solving the system of equations 

2N 

- E 9sC^-'J = 0, n = -M,...,-l,0,l,...,M, (A.2) 

s=-2M 

E c^*ct'^ = h (A.3) 

n=-M 

that results, by linearization [37], from (jS]) and from the truncated normalization 
condition 

E |Cn|^ = 1- 
n=-M 

Obviously, the values of with m < —Ai oi m > J\f entering in some of Eqs. (1A.2I1 
should be taken equal to zero, in accordance with the truncation done. The iteration 
process is stopped when the difference between consecutive solutions is satisfactory. The 
outcome may serve as starting point for a new iteration process with larger values of Ai 
and A", to check the stability of the solution. 

The iteration process just described needs initial values {z/*^°\ c^^-*} not far from the 
true solution. The two different values of u can be obtained from the two eigenvalues 

Xj = exp{2i7iUj) , j = 1, 2 , (A.4) 

of the circuit matrix C [38j for the singular point at ^ = 0. The entries of that matrix 
can be computed by numerically integrating Eq. ([3]) on the unit circle, from z = exp(O) 
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to z = exp(2i7r), for two independent sets of initial values. If we consider two solutions, 
Wa{z) and w^^z), obeying, for instance, the conditions 



then 



Wb{e 

Cu 
C21 



0\ 



1, 
0, 



0, 
1, 



w„ 



C12 
C22 



and 



2i'K 



In 



+ 4Ci2C*21 



(A.5) 



Cu + C22 ± V {.C11—C22) 

The two signs in front of the square root produce two different values for unless the 
parameters Qs in Eq. (|3]) be such that (Cu — C22)^ + 4Ci2C2i = 0, in which case only one 
multiplicative solution appears, any other independent solution containing logarithmic 
terms. The ambiguity in the real part of v due to the multivaluedness of the logarithm 
in the right hand side of (lA.SP reflects the fact already mentioned that the indices v are 
not uniquely defined. Notice that 

Ai A2 = det C = W[wa, Wb] = 1 

and, therefore, 

I'l + 1^2 = (mod 1). 

This may serve as a test for the integration of Eq. (|3]) on the unit circle. 

Although Eq. ( ]A.5I) is exact, the Cmn are obtained numerically and the resulting 
values of u may only be considered as starting values, for the Newton iteration 
process. As starting coefficients c!^^ one may use the solutions of the homogeneous 
system 

2N 

SO) _n 



(^ + ^(o))(^_l + ^{o))c{o)_ ^ g,c^'>, = 0, n = -M,... ,-1,0,1,. ..,Af, (A.6) 

s=-2M 

with the already mentioned truncated normalization condition 



E 



.(0)1 2 



1. 



(A.7) 



-M 



Appendix B. Wronskians of Floquet and Thome solutions 

The connection factors of the Floquet solutions with the Thome ones are given, by 
Eqs. (|2T]l and (|22l) . as quotients of two Wronskians. Those in the denominators can be 
obtained immediately and were given in Eqs. ( I23ll and (IMl) . Direct computation of the 
Wronskians in the numerators must be discarded for the reasons pointed out at the end 
of Sec. 3. The purpose of this Appendix is to give a procedure to compute them. The 
idea is to find, for each one of the needed Wronskians, two functions, one proportional to 
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the other, the proportionahty constant being that Wronskian. Comparison of analogous 
terms in the asymptotic expansions of those functions allows one to obtain the required 
Wronskian. This idea has already been exploited in the solution of the Schrodinger 
equation with a polynomial potential [321 133] . 

Let us consider the Wronskian of one of the Floquet solutions Wj {j = 1, 2), given 
in ([5]), and one of the Thome solutions at infinity Wk {k = 3,4), given in (j6]). We find 
convenient to introduce auxiliary functions 

Mj-fc = exp( — —z jwj, Mfc = exp(--^2 jwk, J = 1,2, fc = 3,4. (B.l) 
Obviously, 

WKfc, Uk] = exp (^-^ W[w„Wk]. (B.2) 

An asymptotic expansion of the left hand side of this equation can be calculated by 
using the definitions fIB.ip and the expansions (E]) and It becomes 



N-l 



where we have denoted 



V 



, = expl t/;,, j = l,2, A; = 3,4, (B.4) 



p 



Sk=Y.(^rn,kZ-'^^^\ fc = 3,4. (B.5) 

m=0 

For the newly introduced function Vj^k-, a convergent Laurent expansion 

oo 

%fc= E Cn,,,fc^"+'^^ (B.6) 



n=— oo 



can be obtained as Floquet solution of the differential equation 



dz"^ \ ^ I dz 



+ 



2N N-l N-l 



E ^.^^+Eb-l)«M^'- E«M^M |^.> = 0. (B.7) 

s=-2M p=2 \p=l 



Then, by using Eqs. (IB.Sp and (IB.6P in (IB.3p . we obtain 

oo 

>VKfc,H~ E 7i'''^ (B.8) 

n=— CO 

where 

oo / N-l 

In'''^ = E '^'TJ.'fc ( '^N,k Cn+m+l-N,j,k + 2 E '^P.^ Cn+m+l-p,j,k 
m=0 \ p=l 

- {n + 2m + l + Uj - /ifc) Cn+m+i,j,k ■ (B.9) 
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The value of the W[wj,Wfc] can be immediately obtained if we are able to write an 
asymptotic expansion of exp (^—ajsi^kZ^/N^, in the right hand side of flB.2l) . with the 
same powers of z as the expansion (IB.SP of the left hand side. With this purpose, we 
construct formal expansions 

of exp (^—aN,kZ^ /N^- Such expansions are but particular forms of the so called 
Heaviside's exponential series 

introduced by Heaviside in the second volume of his Electromagnetic theory (London, 
1899) and which, as proved by Barnes [39], is an asymptotic expansion for arbitrary 
6 and |arg(t)| < tt. Expansions of this kind have been already used by Naundorf 
[27] in his treatment of the connection problem, from which our method is a convenient 
modification. It becomes evident that, for any set of constants {k^''^''} (L = 0, 1, . . . , A^— 
1) satisfying the restriction 

N-1 

^4^»=WK,^d, (B.12) 

L=0 

one has from flB.2l) 

N-1 
L=0 

By choosing for the ^j^''^'' in the expansions S^''^^ the values 

= {iyj+f^k + L)/N, (B.14) 

a comparison, term by term, can be done of the resulting expansion in flRl3|) with that 
in (IB. 81) . One obtains in this way 



, 1 , r(i,fc)x -^nN+L^ (B-15) 

r(n + 1 + ') 



for any positive integer n large enough to satisfy 



2N 



\{n + Uj){n-l + Uj)\> \9s\- 

s=-2M 

By substituting in flB.121) the values of k^'^^ obtained from flB.lSP one has finally 



Tin ^^ -I- 

WK, «^.] = E ^ " ( j) 7i^'|., (B.16) 



^0 (-a^^,/iV)"+^- 

where the minus sign in front of a^^k is to be interpreted as an odd power of e*'^ or e~*'^ 
so as to have | arg(— a^r fc 2;^)| < vr. In the physical problems, one is interested on the 
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connection factors on the positive real semiaxis, that is, on the ray arg(2;) = 0. In this 
case, arg(— ctAT^s z^) = 0; then, Eq. flB.16P gives for one of the numerators in fl2T]) 



^-1 rfr, J- 1 -L A^-^'^h 

L=o {\aN,3\/Nr^ ' 

and the connection factor Tj4 can be obtained immediately from ( 1B.17I) and ( l23l) . For 
the computation of the other connection factor, instead, one has to bear in mind 
the fact that the semiaxis arg(2;) = is a Stokes ray for Tj^^. Actually, for z on 
this ray, | arg(— aAr^4 2;^)| = vr and the expansions ( IB. 101) would not correspond to 
exp (^—ai\[,k z^ /N^. Following the common practice, we define Tj 3 on the Stokes ray by 
the average 

T,,z = \iTl, + T-,) (B.18) 
of its values in the regions separated by that ray. This corresponds to define 

>VK-,H = \ (WK-,W4]+ + WK-,W4]-) , (B.19) 

an average of the Wronskians for z slightly above and below the positive real semiaxis. 
One obtains in this way 

WK-, H = (-1)" E (*''^ ^) ^ 'J ^^unIl , J = 1, 2 , (B.20) 

and the expression of Tj^a follows then from flB.20p and f l23|) . 

The procedure to obtain the Wronskian of one of the Floquet solutions, Wj 
(j = 1,2), and one of the Thome solutions at the origin, wi {I = 5,6) given in (^^, 
is analogous to that just described. The auxiliary functions are now 

Mj-, = exp I \wj, "' = ^^pI~^^ )^'' J = 1,2, / = 5,6.(B.21) 
Then, 

W[u,^u ui] = exp (^-^ z-''^ WK-, wi]. (B.22) 

By using the definitions (IB. 211) and the expansions ([5]) and ([7]) one obtains an asymptotic 
expansion of the left hand side, namely, 

>VK,,, ui] ~ (^-Pm,i z-''-' - 2 e' z"^-'^ v,,i - ^ j Si + v,,i ^, (B.23) 
where we have denoted 

(m-1 n \ 

v,,i = e^^r£^^'j'^'' J = 1,2, / = 5,6, (B.24) 

oo 

Si=Y.h^^iz^+p^ , 1 = 5,6. (B.25) 

m=0 
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oo 
n=— oo 

for Vj^i can be obtained as Floquet solution of the differential equation 
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(B.26) 



1^ 



' M-l 



+ 2z \ 



-9 



'hi 



q=l 
( 2N 



+ 



dz 

M-l 



'M-l 



s=-2M 



q=l 



9=1 



From flR23|) . by using Eqs. flR25|l and flR26|) . we obtain 



(B.28) 



where 



in 



oo 

E 

m=0 



Af-1 



/3m,/ Cn-m+l+M,j,l ~ 2 ^ Cn-m+l+q,j,l 

9=1 

+ (—77, + 2m — 1 — Uj + pi) c 



■n—m+l,j,l 



(B.29) 



In order to obtain an asymptotic expansion of exp \^—l3M,iz~^^/Mj, in the right hand 
side of ( 1B.22I) . with the same powers of z as the expansion ( IB. 281) of the left hand side, 
we construct M formal expansions analogous to those in fIB.lOl) . 



oo ( 



E 



(3,1) 



Let us now consider M constants {/t^''''} [L = 0, 1, 

M-l 

Then, 

M-l 
L=0 

If we choose for the (jj^'''' in (IB.30P the values 
6^1'^ = i-u,-pi + L)/M, 



L = 0,1,..., M-l. (B.30) 
M-l) such that 

(B.31) 
(B.32) 



(B.33) 

comparison of the resulting expansion in (1B.32P with that in (1B.28P allows one to write 



that, substituted in (IB.3ip . gives 



W[w„wi] = J2 



L=o {-/3m,i/M) 



n+5\ 



l-nM-Li 



(B.34) 



(B.35) 
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where the minus sign in front of should be replaced by an odd power of e*'^ or e ^'^ 
such that I aig{—PM,i < ^r. For z on the positive real semiaxis, one has 

Fl'-n -I- 1 -I- A'^-^'^h 

^=0 (|/3m,5|/M)"-*-^- 
and, analogously to Eq. (1B.20I) . 

The connection factors Tj 5 and Tj g are then immediately obtained from Eqs. ( |22|) by 
using (EID, f lR36l) and f lB37l) . 
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